CoUisionless Absorption of Intense Laser Beams by Anharmonic Resonance 
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Two decades after the invention of chirped pulse amplification the physical mechanism of collisionless ab- 
sorption of intense laser radiation in overdense matter is still not sufficiently well understood. We show that 
anharmonic resonance in the self-generated plasma potential of the single plasma layers (cold plasma model) 
or of the individual electrons (warm plasma), respectively, constitutes the leading physical mechanism of colli- 
sionless absorption in an overdense plasma. Analogously to collisions, resonance provides for the finite phase 
shift of the free electron current relative to the driving laser field which is compulsory for energy transfer from 
the laser beam to any medium. An efficient new scenario of wave breaking is also indicated. 
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PACS numbers: 52.25.0s, 52.38-r, 52.38.Dx, 52.35.Fp 

Collisionless absorption of intense sub-ps laser pulses in 
solid targets up to 70% and beyond is possible [1]. The results 
found their confirmation in numerous particle-in-cell (PIC) 
||2,[3t] and Vlasov simulations P, Since the beginning in- 
tense search for physical mechanism supporting such high ab- 
sorption values in the absence of collisions started and a whole 
variety of effects was made responsible for the observations. 
Pioneering work in this context has been done by several re- 
searchers. W. KruerandK. Eastabrook introduced the concept 
of j X S-heating in 1985 [6] . They realized for the first time 
the relevance of electron motion in laser beam direction for 
heating. A further, very significant step forward towards un- 
derstanding strong collisionless absorption was made in 1987 
due to F. Brunei by recognizing the role of the collective self- 
generated electric field in the plasma under oblique incidence 
or the action of the j x S-forrce, respectively [7]. Finally, P. 
Gibbon discovered that a certain fraction of electrons pushed 
out into the vacuum do not return to the target with the rhythm 
of the laser oscillations. For this phenomenon, because of in- 
dicating irreversibility aspects, he coined the term "vacuum 
heating" in 1994 which subsequently has become very pop- 
ular |2, 's^. In our opinion it was unfortunate that this phe- 
nomenon did not find the systematic attention of researchers 
it deserved. So far also quantitative analytic absorption mod- 
els have been elaborated, like anomalous skin layer absorption 
Jgt], stochastic heating ifioll . or Landau damping |11], which 
however yield absorption values approximately an order of 
magnitude lower than the values above we are interested in 
in this context here. For a more detailed analysis see il2ll . 
According to some researchers' view one of the phenomena 
makes the contribution to collisionless absorption, according 
to others althogether contribute. In summary, collisionless ab- 
sorption in overdense matter is well confirmed by experiments 
and simulations but not understood physically to the desired 
extent. 



We introduce the new concept of anharmonic resonance and 
show that it is the leading physical mechanism responsible for 
collisionless absorption. In a first stage we give a descrip- 
tion of the idea for the hurried reader; without loss of physi- 
cal insight he may skip all formulas. Then we show in detail 
how anharmonic resonance works with a single plasma layer 
Subsequently we study numerically the interaction of multi- 
ple layers coupled together and, with the aid of test particle 
simulations, we present the proof of our assertion for realis- 
tic interaction, with the limitation to plane geometry for being 
basic in many respects. 

Formally, the search for collisionless absorption reduces to 
the question where the phase shift between the monochro- 
matic electric field of the laser E ~ cos Lot of frequency lu 
and the electron current density j ^ sin(wi + (j)) in Poynt- 
ing's theorem comes from in the absence of collisions since, 
with (j) = 0, j ■ E sin ojt cos uut = 0. The most efficient 
collisionless mechanism producing (f> ^ Ois resonance. How- 
ever, resonance absorption atuip^uj in the plasma has been 
categorically excluded so far by the scientific community. The 
insight that this statement is incorrect, because limited to the 
harmonic oscillator only, is the key to the solution. In fact, 
when the electrons of a plane plasma layer of thickness d are 
displaced by an amount ^ = d or larger their oscillation eigen- 
period Tq lengthens from 2Tr /uip to larger values and at a suf- 
ficiently big excursion ^ becomes equal to the laser period 
T = 27r/aj, i.e., it enters into resonance. This happens already 
at low laser intensities. The transition from non-resonant to 
resonant state is irreversible and is accompanied by the de- 
sired well-known phase shift of \(f>\ = tt. This behavior per- 
sists also in the real plasma with crossing of multiple oscil- 
lators (breaking) and mutually influencing each other thereby 
lowering the single layer resonance threshold. Owing to the 
non-integrable (chaotic) nature of the governing equations al- 
ready in the simplest case this final part of the assertion can 
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be proved only with the help of numerical methods. 

Now we proceed to the quantitative proof. A plane, fully 
ionized target is assumed to fill the half space a; > 0. A plane 
wave E{x,t) — EQexp[i{k ■ x — ojt)] in p-polarization in 
y-direction, wave vector k and frequency uj, is incident from 
— cx) under an angle a onto the plasma surface (Fig.[T^). Af- 
ter applying a Lorentz boost vq — c sin a the wave impinges 
normally. The target is thought to be cut into a sufficient 
number of thin parallel layers, each of which is exposed to 
a driving force F acting along x of magnitude cvqE and fre- 
quency Lu, and an additional component originating from the 
Lorentz force of the oscillatory motion along y of frequency 
2uj (e electron charge, B magnetic field of the laser). With 
the electron displacement ^ in x-direction and immobile ions 
of density uq (corresponding to ojpo) the motion of a single 
layer (Fig.[T}5) in the non-relativistic limit is well determined 
by the potential y = me{ujpod/2)'^{l + (^)^^^, ( = 
electron mass (derivation is elementary). Hence, the plasma 
oscillator is governed by the dimensionless equation 

where r = Wpo^ and D = 2F/{Lo^^d). The eigenperiod 

To = y 4[2(£ - F(e))/me]-i/2 (2) 

for a fixed energy S, at high excitation and ;» is 
given by Tq ~ 8(ijjpod)~^/^^Q^^. The eigenfrequency cjq = 

2-1: /To = (7r/4)(tjpQ(i)^/^^Q decreases with increasing os- 
cillation amplitude in contrast to the linear harmonic oscil- 
lator To make the essence clear we concentrate on large an- 
gles a and set D — ao{t) cos wt, ao(t) amplitude. At = ^ 
the single plasma layer is resonantly excited at amplitude 

= — {TTU!pQ/Au!)'^d and t = tr- Under the assumption 
that the oscillator is driven into resonance during half a laser 
cycle, i.e., when the driver amplitude Eg = meUJpQd/ [Ae], 
with d — 0.1-0.2 nm and UpQ — 2y. 10^^ s~^ this happens at 
the laser intensity / = 10^^ Wcm^^. 

Under a weak driver D{t) in ([T]l the eigenfrequency tJo is 
much higher than lo and the excursion C(t) ~ cos{Lot + ip) 
follows adiabatically D{t), i.e., 1^9 = 0. In other words, C,{t) 
moves in phase with D{t). At large displacements the restor- 
ing force weakens and loq <^ lo, thus imposing \lp\ — tt, like a 
free electron oscillating in the laser field. At an intermediate 
C-value the phase shift must have taken place with the conse- 
quence that the product (D transits from ^ sin r x cos t to 
cos^ r. It is the resonance region. Only there, in the so-called 
neighborhood of stationary phase |13] resonance occurs and 
irreversible energy absorption by the oscillator is possible. A 
necessary condition for stationary phase is the existence of a 
point with \(p\ = 7r/2. Slow transitions through resonance can 
be visualized geometrically with the help of the Cornu spi- 
ral as explicitly shown for uJo{t) in another context in lll4ll . 
In order to get rid of such WKB-like excitation Eq. ([T]i is 
solved numerically with a smooth A^-cycle sin^-shaped driver 




FIG. 1 : (a) Oblique incidence of parallel laser beam of wave vector 
k in the lab frame and k' || x in the system boosted by i^o = csin a. 
The overdense target is cut into layers of thickness d. (b) Large elec- 
tron displacement ^ in plasma layer (LHS) and anharmonic closed 
(dotted) and semi-open (dashed) potentials. In wider potential than 
harmonic (solid) the frequency decreases with increasing ^. 
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FIG. 2: Excitation e{C) = /2 + V{C) of the oscillator from 
Eq. l[T| by the driver D{r) = Do sm'^[ujT/{2Nujpo)] coslot/lopo, 
uj/iOpo = 0.1, = 20, (a) below resonance (Do = 0.921), (b) 
above resonance (Do — 0.923) with an 43 times higher final en- 
ergy gain Ef than in (a) although the driver strength is changed only 
by 0.2%. The potential and the resonant energy level are indicated 
dashed and dotted, respectively. Note the different scales in a, b. 

D = ao (i) cos tji = D o sin^ [ujt/ {2 N)] cos ujt. In Fig. |2^ 
with a driver Dq = 0.921 the layer remains below resonance, 
the energy gain ef — e{2NTTU!pQ/uj) after the pulse is neg- 
ligible. By an increase of the driver of only ADq = 0.002 
resonance takes place. Due to dephasing above resonance the 
energy content shows the typical periodic variations attenuat- 
ing in time towards a finite value when the pulse is over. Much 
energy is stored now in the oscillator (see the horizontal orbits 
in Fig.|2]3). Compared to case (a) the energy gain is increased 
by a factor of 43. 

Figure[3]is of particular relevance. It shows the driving field 
D{t), the displacement C,{t), and the energy gained e(r). At 
position 1 the oscillator is entering resonance (e starts increas- 
ing, cjQ > lo), D and C, are in phase; at 2 it is leaving resonance 
{e starts decreasing, o^o < lo), D and ( are dephased by tt. Po- 
sitions I and II (points of stationary phase) indicate maximum 
energy gain and maximum energy loss, D and ( are dephased 
by ±7r/2. Thus, the resonance signature is preserved in a rapid 
transition. The phenomenon repeats in the second maximum 
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(non-integrable, chaotic), yet elementary Hamiltonian 
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FIG. 3: Driver D{t), excursion ((t) and absorbed energy e(r) vs 
time (in driver cycles) for the case in Fig. |2] Each time resonance 
is crossed undergoes a phase shift by n. Bold lines I, II indicate 
instants of phase lag are ±7r/2 between D and , i.e., maximum 
energy gain and loss (points of stationary phase); modulations in e 
originate from the oj + a;o spectral component. 
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FIG. 4: Resonance dynamics of layer #32 as a function of laser cy- 
cles. Position kx (a), k wave number, absorbed energy e (b) and 
phase ip between velocity v and driving laser field E — Eo{x) cos uit 
(mapped into the interval [— vr, 0]) (c); e in units of ponderomotive 
potential Up. Passages through — 7r/2 are indicated by I, II, and 
1 and 2 (see text). Eo{x) is calculated selfconsistently and corre- 
sponds to laser intensity I = 3.5 x 10^** Wcm^^. After resonance 
at 3.5 laser cycles the layer is pushed back into the target with high 
velocity ("disruption"). After crossing the opposite target surface the 
layer is substituted by a new layer Note E ^ —D. 



of e, etc. The definition of resonance as the points of station- 
ary phases is the natural extension of the concept of resonance 
to anharmonic and nonlinear oscillators. An equivalent crite- 
rion to be used later is this: (i) the half widths of the local 
maxima of e(r) is not much shorter than half a driver cycle 
T = 27r/aj (in Fig.|2]it is 1.5T), in contrast to a collisional 
event which is almost instantaneous, and (ii) e(r oo) > 0. 

The motion of N layers of free electrons and ions is de- 
scribed for arbitrarily large oscillations by the non-separable 



fe=l y k'^k 1=1 



Vki - D{T)Ck 



(3) 



withpfe = dCfc/dr, Vkk' = [1 + (Cfe - Cfe')']'/', Vki = 
-[1 + iCk - Co/)']'/', a - '2xk/d, Coi - 2ai/d. When one 
of the layers is driven into resonance it starts moving opposite 
to the coherently moving non-resonant layers thereby crossing 
one or several adjacent oscillators. This is a new scenario of 
very effective wave breaking not described in the literature so 
far. We give it the name of resonant (wave) breaking. It means 
loss of coherence and leads to flattening of the collective po- 
tential and, in concomitance, to a reduction of the resonance 
threshold. As a representative case we study the dynamics of 
a 100 times overdense target, subdivided into 120 layers of 
d = 0.125 nm each, on which / 3.5 x 10^^ Wcm^^at 
A = 800 nm is impinging. The typical scenario is as follows: 
After being pulled out into the vacuum and oscillating there 
for some time the layers are pushed back in a disruption-like 
manner into the target (formation of jets). When the layers 
leave from the back of the target they are replaced by new lay- 
ers with zero momentum (cold return current). First indica- 
tion of resonance: After 100 laser cycles more than half of the 
now more than 2700 layers have gained energies exceeding 
their quiver energy Up. The power spectrum shows a plateau 
between lUp and the cut off at 6C/p. To show the occurrence 
of resonance explicitly the phase of each layer with respect 
to the driving laser field was investigated. A typical example 
is shown in Fig. |4] with layer #32, LHS trajectory and driver 
field (a), middle absorbed energy (b), RHS phase ip of veloc- 
ity V ^ sm{ojt + (fi) with respect to the driver E ^ cos ust 
(c). Over T/2 there is a continuous and smooth transition 
of (p through — 7r/2att/T = 2.8 (I) with a simultaneous 
strong increase in the absorbed energy (b) and the excursion 
X, with following disruption of the layer at t/T — 3.6. The 
change of (p is clearly seen also in (a). Another resonance 
of the same kind is found at t/T = 8.8 (II). Other two pas- 
sages of (f through — 7r/2 at t/T — 5.7 (1) and 8.1 (2) show 
rapid fluctuations and hence do show almost no energy gain 
and no disruption [see (b), (a)]. Transitions of this latter kind 
are morphologically clearly distinguishable from the former 
case, and for none of the 120 layers they are able to acceler- 
ate them across the target. This proofs in an impressive way 
that our definition of resonance as points of stationary phase 
and the resonance criterion (i)H-(ii) are correct and meaning- 
ful. We conclude that in the cold plasma model all absorption 
is by anharmonic resonance. The layers disrupt in the chaotic 
order6, 5, 4, 3,2, 15, 17, 16, 9, 11, 10, 8, 12, 24,31,34, 28, 
32, 14, 30, 1, 33, 45, 26, 36, 27, etc.; layer 113 disrupts before 
front layer 0. This is one of the fundamental differences in 
the dynamics in comparison to [16]. More essential is that the 
Hamiltonian Eq. (3) acts on the single layer like a half open 
potential, as sketched by the dashed line in Fig. [T]^. We tested 
this explicitly and found that the degree of absorption does not 
depend much on the hight of the potential barrier, even when 



4 




x/ urn 

FIG. 5: 1D3V PIC simulation of a Nd laser beam interacting with 
a thick overdense target with mobile ions in the boosted frame (pa- 
rameters see text). Regular shadow structure: laser field, LHS black 
trajectory: displacement x(t) of test electron #4, RHS black trajec- 
tory starting from OrrieC: corresponding momentum px, white line: 
total electric field at position x{t). Resonance (strong momentum 
increase and phase shift and disruption are impressive. No 

disruption without preceding resonance ever seen in any test. All 
resonances look very similar to each other 



this latter lies considerably below the resonance level of the 
closed potential. Only when it is set equal to zero as in |7] 
the absorption almost vanishes because of the absence of res- 
onance. The half open potential with finite threshold towards 
the target interior is clearly seen in our PIC simulations with 
the PSC code. 

The final step of the proof , i.e., high degree of absorption 
as mentioned in the introduction jHUIHB] ™d absorption it- 
self being accomplished by anharmonic resonance, at present 
can only be based on computer simulations. For this purpose 
a 10^ particles PIC simulation with the PSC code OlSIl in its 
collisionless mode under 45° iiTadiation is performed. In the 
boosted frame a Gaussian Nd laser beam of / = 10^^ Wcm^^ 
and halfwidth of 26 fs acts on a plane 80 times overcritical 
7.3fim thick target. The orbits of 200 test electrons equally 
distributed over the target are followed during 15 laser cycles 
of = 5 fs. Again, the outcome of the statistics is over- 
whelming: all test electrons interacting with the laser field 
are resonantly accelerated and disrupt nearly immediately. In 
Fig.|5]the time history of test electron is depicted. The electron 
enters the laser field (shadowed interference pattern) and inter- 
acts resonantly (see momentum) and escapes into the target an 
instant later with 5 times Up. All resonant orbits look nearly 
identical to each other. The resonance character is ensured 
by the high energy gain and the phase shift in comparison to 
the total laser field (white line superposed). The shadowed 
fine structure right of the laser field is due to plasmons. The 
simulation tells also important details on the heating mecha- 



nism. The primary effect is the generation of the fast spec- 
trum by resonance. They generate "solitary", i.e. non-Bohm- 
Gross plasmons in the dense interior which, in turn, heat cold 
electrons by a mechanism resembling Landau damping. This 
completes our last part of the proof. It becomes clear now 
why perturbative theories when starting in zero approxima- 
tion from straight orbits, e.g., anomalous skin effect, fail to 
explain strong absorption because anharmonic resonance is 
outside the validity of a linearized treatment of standard type. 
We want to point out that this mechanism is active also in long 
fs or ps pulses when profile steepening is so strong that no lin- 
ear resonance can take place. 

In summary we have addressed the leading physical mech- 
anism of collisionless absorption and, in particular, we have 
discovered the phase shift between driver field and electron 
current, indispensable for absorption, as a resonance effect. 
The discovery may be viewed as the extension of the well- 
known fact that a single point charge cannot absorb a photon 
unless it resonates in an outer potential. As a byproduct we 
have found a new, very efficient scenario leading to (wave) 
breaking. Our results will have a major impact on the fur- 
ther progress in the theory of laser-dense matter interaction. 
So, for example, on the basis of anharmonic resonance the ap- 
pearance of a hot temperature lying higher, often considerably 
higher, than the mean electron quiver energy and the fact that 
the Maxwellian-like tail of their energy spectrum is filled up 
nearly instantaneously (fs time scale) finds a very natural ex- 
planation. The latter phenomenon, never decently discussed 
in the literature, is particularly surprising as the fast electrons 
do almost not interact together Finally, as the overwhelming 
majority of oscillatory motions in nature are anharmonic, the 
harmonic oscillator being the great exception, the model de- 
veloped here will find its application in various other fields of 
fundamental and applied science , e.g., formation of cracks by 
fatigue under oscillatory stress, and in catastrophe theory. In 
the specific field of high power - overdense matter interaction 
the main relevance of our finding we see in the possibility to 
tailor the electron spectrum for various applications (electron 
and ion acceleration, fast ignition, etc.) by designing targets 
properly, for instance by choosing carefully their thickness. 
The energy spectrum of thin targets is more energetic owing 
to multiple resonances than the thick target spectrum. Finally, 
the simulations have also revealed that in the latter case the 
spectrum is subject to continuous metamorphosis in time, an 
aspect which may play an important role in fast ignition. 
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